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Abstract 

Effective equations often provide powerful tools to develop a systematic under- 
standing of detailed properties of a quantum system. This is especially helpful in 
quantum cosmology where several conceptual and technical difficulties associated 
with the full quantum equations can be avoided in this way. Here, effective equations 
for Wheeler-DeWitt and loop quantizations of spatially flat, isotropic cosmological 
models sourced by a massive or interacting scalar are derived and studied. The re- 
sulting systems are remarkably different from that given for a free, massless scalar. 
This has implications for the coherence of evolving states and the realization of a 
bounce in loop quantum cosmology. 

1 Introduction 

While initial stages of the expanding branch of our universe are becoming better and better 
understood, its potential history before the big bang remains open to many debates. By 
now, several scenarios have been suggested which, building on quantum gravity effects 
rather than purely classical gravity, could give rise to a non-singular evolution at all times 
[H El m m [HI E] . The era preceding the classical big bang singularity, as well as the quantum 
transition through it, can then be relevant to explain observed features of our universe and 
must therefore be thoroughly understood. Most of these scenarios realize the simplest 
possible non-singular behavior where an isotropic collapsing universe reaches a smallest, 
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non-zero volume before it bounces back into an expanding branch. In addition to isotropy, 
one often assumes a special form of matter, such as a free massless scalar making quantum 
cosmo logical models solvable. Whether or not this results in a robust, reliable scenario 
for a complex universe can only be seen by embedding such simple models in a more 
general class, at least perturbatively, including matter interactions or inhomogeneities. 
From a conceptual viewpoint, this provides a stability analysis of the original model; from 
a phenomenological viewpoint it will allow one to construct realistic scenarios with all 
ingredients relevant for observations. 

Loop quantum cosmology gives rise to a potentially general mechanism which removes 
space-time singularities [3 El E] • In this general form, not much is said about the precise 
regime around the big bang and the state before, and especially at, the classical singularity 
might have been highly quantum. This, in fact, agrees with general expectations. It 
thus came as quite a surprise when detailed properties derived for an isotropic model 
in loop quantum cosmology sourced by a free, massless scalar indicated that the state 
remained semiclassical to a high degree throughout its evolution, before and after the big 
bang [To]. This model was analyzed because the absence of a potential allows a direct 
calculation of observables and an implementation of the physical inner product. This, in 
turn, made it possible to evaluate numerical solutions to the underlying evolution equations 
for physical information [TTl 112] . When such special mathematical simplifications, beyond 
those already implied by assuming space-time symmetries, are used it is never clear whether 
detailed properties are general or specific to the chosen model. This especially applies 
to the semiclassical aspects which were unexpected from the point of view of general 
quantum systems whose states usually spread. An analytical analysis of a related and 
exactly solvable mode0 [TU [15] clarified this issue and showed that the semiclassicality 
properties were indeed very special: the bounce is described by a solvable model in which 
expectation values and fluctuations (or higher moments of a wave function) do not couple 
to each other. Thus, while the state evolves, its expectation values are not affected by its 
spreading, and the spreading is completely independent of other parameters of the state. In 
other words, the system is analogous to a harmonic oscillator in quantum mechanics whose 
states can remain coherent forever while following the classical trajectory. As quantum 
systems go, this high degree of solvability is a rare property well-known from the harmonic 
oscillator or free quantum field theories but not corresponding to a general feature. 

In this model, expectation values, spread parameters and higher moments of a state 
evolve independently of each other and can be determined exactly. Thus, while the wave 
function evolves and possibly spreads and deforms, this does not influence the trajectory 
followed by its expectation values. But when details of the model are changed, solvability 
can no longer be maintained and complicated coupling terms between expectation values, 

^Solvability requires one to drop some of the ingredients of a loop quantized Hamiltonian which, how- 
ever, are not relevant for states of a universe with large matter content. The solvable model is thus not 
exactly the same as a loop quantization or that used in the numerical studies of [TT] but still allows one to 
analyze the bounce. Note that also in [TT] the dynamical equation was adapted for the numerical purposes 
and differs from what one would obtain in a loop quantization |13j . None of these changes matter for 
properties of the bounce of a large universe. 
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spreads and higher moments arise. This leads to quantum corrections to the classical equa- 
tions which can often be captured in effective equations [HI [T71 [18]. In this paper, we set 
out to study such effective equations for an isotropic model in loop quantum gravity with 
a massive and self-interacting scalar in order to test the genericness of properties of the 
exactly solvable model. While the solvability of the exact bounce model certainly makes 
its properties very special, it also provides the starting point for a systematic perturbation 
analysis. This is analogous to free field theories used as the zeroth order in perturbative 
quantum field theory. Adding a potential as done here is the simplest generalization in 
that it does not introduce additional degrees of freedom. (Initial steps to include inho- 
mogeneities have been done in [19].) However, it also leads to difficulties since the scalar, 
instead of being a monotonic function of coordinate time as in the free case, has itself non- 
trivial dynamics. As discussed in more detail below, a monotonic scalar is used crucially 
in the solvable model since evolution of the universe volume is described relationally with 
respect to the scalar. If the scalar is no longer monotonic, it cannot be used as a global time 
coordinate, meaning that the model with a potential can only be used for a finite range of 
time. Nevertheless, we will see that this allows one to test whether quantum back-reaction 
effects remain negligible or whether they can have a strong effect on the behavior of a 
bounce. 

It is often assumed that a potential term of a scalar field cannot have a sizable effect near 
an isotropic classical singularity because its classical form a^V{(j)) tends to zero for a — >■ 
while the kinetic term ^a'^p"^ in canonical variables has a diverging pre-factor. The general 
behavior certainly depends also on how and behave when the singularity is approached, 
but if the potential is small does not change much, being a constant of motion in 
the free case. Thus, the kinetic term is expected to dominate and a non-zero potential 
should not change much of the bounce observed in free models. However, this reasoning 
overlooks the quantum nature of the problem. It is based on a classical Hamiltonian 
and the evolution equations it implies. If quantum gravity is used to argue for the free 
bounce, quantum aspects must also be taken into account for interactions implied by the 
potential. Here, details of the free bounce related to solvability become important. The 
solvable bounce model demonstrates that a state which starts out semiclassically at large 
volume and evolves to smaller volume stays semiclassical and then, in this semiclassical 
form, bounces when quantum geometry effects set in. The model also shows that there 
do exist states, though not semiclassical ones, which do not bounce [15]. This does not 
alter the conclusions about singularity removal in the free model since one always has 
to use the boundary condition that a state is semiclassical at large volume (at least at 
one side of the bounce). Under this condition, every solution of the free model bounces. 
But this is no longer guaranteed with interacting matter. Any realistic statement about 
singularity removal has to involve evolution over long stretches of time: a universe which 
is semiclassical at large volume, like our own, must be shown to remain non-singular after 
a long time of backward (or forward) evolution. This does not play a large role for a 
solvable model since there are no quantum back-reaction effects of an evolving state on its 
expectation values. But a non-zero potential introduces such back-reaction effects which 
must be considered even if they are small for a perturbative potential. A reliable statement 
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about non-singular behavior must involve a proof that a state starting out semiclassically 
at large volume does not spread and deform too much so as to avoid a bounce, which may 
well be possible since the free model allows non-semiclassical states which do not bounce. 
In general, the bounce regime, or whatever regime is valid around the classical big bang 
singularity, is expected to be of a highly quantum nature; arguments based solely on states 
which are semiclassical in this regime, by choice or by the lack of quantum back-reaction 
in a specific model, cannot be complete. 

In this paper, we provide effective equations taking into account quantum back-reaction 
in the presence of a perturbative potential and provide an initial investigation. This is a 
first step of a stability analysis of bounce models in loop quantum cosmology. We will find 
several new properties and arrive at a set of effective equations of a new type requiring 
the presence of independent quantum degrees of freedom. Even though by construction we 
need to require the potential to be small and a perturbation to the kinetic term, several 
new qualitative features can be seen. The analytical, though approximate, analysis allows 
more general conclusions than numerical investigations, which are often prone to hidden 
assumptions, would show. 

2 Description of the models 

We start with the Friedmann equation in an isotropic and homogeneous background with 
a scalar field 0. It can be written as a constraint C = with 

in terms of canonical variables \p\ = a^, c = da/dr (using proper time r) for the metric 
and = a^d(f)/dT as the canonical momentum conjugate to 0. We will be ignoring factors 
of SvrG/S in this article to keep subsequent formulas simple. A convenient understanding 
of dimensions then is to use dimensionless c and while p and 1/V have the dimension 
of length squared. Also the Planck constant h has the dimension of length squared if G is 
ignored in hG = ip with the Planck length £p . The variable p is derived from an isotropic 
densitized triad [20] and can thus in general assume both signs to take into account the 
triad orientation. In this paper, however, it will be sufficient to use only positive values of 
p since our aim is to analyze bounces; a change of sign in p would imply that no bounce 
occurs. (Even if quantum effects imply regular evolution through the classical singularity 
at p = 0, this would not be considered a bounce which conventionally requires a non- 
zero lower bound to volume. "Bounce" is usually also reserved for situations in which a 
semiclassical description remains valid at all times and the universe does not enter a deep 
quantum regime. This is realized in all examples for bounces described so far.) 
The Friedmann equation describes the dynamics of an isotropic space-time 

ds^ = -dr^ + a(r)2(dx2 + dy^ + dz^) 

with the scale factor a, sourced by a homogeneous scalar matter field with potential V{(j)). 
Since ([1]) is a constraint, it generates equations of motion by d//dt = {f,NC} for any 
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phase space function /, referring to a derivative by a time coordinate t which corresponds 
to the chosen lapse function A^. For proper time r, as in the above hne element, we simply 
have = 1. Then, 



dp/dr = (2) 

d0/dr = \p\-'^% (4) 
dp^/dr = -\p\'/'V\<P). (5) 



The first and third equation determine the relation between momenta and time derivatives 
of coordinates, while the other two then result in second order evolution equations for p 
and in proper time. 

2.1 Internal time 

Instead of using a coordinate time r, we can solve ([1]) for p^, 

-^W\ = \pWc'-\p\Vi<l>) =:H (6) 

and view the right hand side as a Hamiltonian generating relational evolution in 0. Rather 
than solving for p{t), c(r) and 0(r) where, due to general covariance the choice of time 
coordinate r does not matter, we thus eliminate coordinate time altogether and look for 
solutions p{(f)) and c(0). We will deal here mainly with a perturbative scheme considering 
the potential to be small in comparison with the gravitational contribution. Thus, we will 
often refer to the expanded Hamiltonian 

H=\p\c-^V{<P) + 0{V') (7) 

up to second order in V. 

Relational solutions in terms of can easily be found explicitly in the potential free 
case since p^ is a constant and thus is a monotonic function of r. Any point in r is 
mapped one-to-one to a value of 0. However, if there is a non-trivial potential such as a 
mass term |m^0^, in general will no longer be monotonic but shows oscillations around 
minima of the potential. Only stretches between two turning points of can then be de- 
scribed in an internal time form, but not all of the evolution that would be accessible in 
r. Moreover, the system is more complicated due to an explicit "time" dependence from 
V{(f)) appearing in H. As we will see, upon quantization this still allows us to investigate 
the effects of spreading and deforming wave packets on the expectation values of a corre- 
sponding quantum system. In particular, interesting phases in early universe cosmology 
are inflationary which, if sourced by an inflaton, typically provides a slow-roll phase where 
the scalar rolls down a fiat potential. In fact, we will see shortly that the potential we deal 
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with perturbatively must be flat so as to provide this setting automatically. This implies a 
long monotonia phase for during which back-reaction effects could build up and become 
important. In backward evolution, this allows us to see if quantum effects could possibly 
prevent or dramatically change a bounce. 



2.2 Wheeler— De Witt quantization 

A further difficulty related to the potential arises at the quantum level. To see this we first 
recall from p4j how the potential-free case is treated quantum theoretically. Instead of a 
time dependent Hamiltonian, this case would give us simply -^b</.| = bc| which, when 
quantized, results in a Schrodinger equation 

Solutions to this equation automatically solve the equation 

-y^^ = p'c2^ (9) 

which (in a certain factor ordering) quantizes ([1]) before taking the square root. Since the 
Friedmann equation, from the canonical point of view, is the Hamiltonian constraint, it 
is indeed quantized first without solving for in canonical quantum gravity. Although 
([9]) is a second order differential equation and thus has more solutions than the first order 
equation ([8]), one usually allows superpositions only of solutions with a fixed sign of energy 
P0 which is chosen positive in taking the root above. Thus, the solution space to iQ 
relevant for the quantum theory should be considered as consisting of two sectors, one of 
which is given by solutions to ([H]) and the other by solutions with the opposite sign of the 
square root. 

A solvable model is then obtained by using the Hamiltonian H = ^{cp + pc) for 
0-evolution which is quadratic in canonical coordinates, just as the harmonic oscillator 
Hamiltonian. This is responsible for the solvability, a fact which, more surprisingly, even 
extends to the loop quantization [H] as recalled in the next subsection. Solvability implies 
that the equations of motion (written in internal time whose derivative is denoted by a 
dot) 

.(a> , ,10) 



ih ih 
form a closed system and can thus be solved easily without knowing the complicated ways 
a state may spread and deform. In general quantum systems, by contrast, such equations 
would form a system of infinitely many coupled equations also involving the quantum 
vSjI"! OjIdIcs 

G»'" = (((C-(C)) — (p-(p))'^)weyl> (11) 

for 2 < n G M and < a < n [Ml [I7j, the subscript "Weyl" denoting fully symmetric 
ordering. Strictly speaking, this is even the case in the system considered here since we 
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dropped the absolute value of the Hamiltonian to make it quadratic. As discussed in de- 
tail in [15], the quadratic constraint still allows one to describe states correctly which are 
sharply peaked at large values of energy if or p^, a case of primary interest for isotropic cos- 
mological models. For the quadratic Hamiltonian, explicit solutions to (|TOl) for expectation 
values c = (c) and p = (p) as well as fluctuations are given by 

cv=o{<P) = Cie'^ (12) 

py=o(0) = C^e-^ (13) 

Gt?=o(0) = ^se'^ (14) 

G^^=o(</') = (15) 

G'^Uct>) = C,e-'^ (16) 

which we distinguish from the perturbative solutions with non-vanishing potential deter- 
mined later by the subscript. Since to the orders analyzed here we are mainly interested 
in fluctuations, we will use only quantum variables at n = 2, and call them G^*^ := G^'^, 
G^P ■= and G^^ := G^'^ for better clarity. 

In the presence of a potential we would, following the same procedure, first quantize to 
obtain a second order equation 

- Y^^= (^^-|p|V(</)))^. (17) 

Now, however, [p^, H]^ Q for a quantization of 

H=\pWc^-\p\V{cP). (18) 

Thus, solutions to 

- -^P^i, = Hi; (19) 
do not solve the second order equation f[T7|) but rather 

^pJV = P4>H^I^ = + [p^, H]^P = H^i/j + [P^, H]^Ij . (20) 

There is thus no strict analog of the first order equation, complicating any analysis. 

However, for the solutions we are interested in here, which includes perturbations 
around solutions of the exact free scalar model obtained for large H, we can ignore the 
commutator term. In fact, due to the square root in H we have, up to factor ordering in 
a precise quantization, 

[p^, H] ~ th — 

ZH 

whose expectation value is small compared to that of H^, also appearing in (120|) for states 
with large H and a not too steep potential. This is exactly the regime we intend to probe 
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by our perturbation theory. Such a regime easily arises without fine-tuning and can thus 
provide general insights: The two conditions |p|^V^(0) ^ p^c^ ~ for the potential term 
to be perturbative and hp^V'{(j))/H <^ for the commutator in (|20l) to be negligible 
require V <^ H^/p^ and hV <^ H^/p^. In the loop quantized model to be discussed next, 
a semiclassical bounce in the free model occurs for H ^ h for the universe not to enter the 
deep quantum regime, and we have p > H . Thus, around such a bounce the conditions 
on the potential are hV <^ 1 and hV <^ 1, such that both the potential and its first 
derivative must be small compared to the Planck mass squared. For a sufficiently small 
and fiat potential, long evolution times can be considered. For the purpose of describing 
those regimes, the first order equation ( |T9l) is a good approximation, and due to large p^ 
we can again drop absolute values on H. 

2.3 Loop quantization 

It is more interesting to study perturbations around the solvable model determined by a 
loop Hamiltonian following again [T3] for the zeroth order solutions. In a loop quantization 
there is no operator for c [21] , and one can thus not use a Hamiltonian of the same form as 
before. Instead, all exponentials e^°"^ with a G M are represented by well-defined operators 
and appear in the Hamiltonian such that the classical expression cp is reproduced in the 
classical, small-curvature limit. This is accomplished by the free Hamiltonian psinc as it 
follows from the loop quantization (up to quantization choices which we can ignore here). 
This Hamiltonian is no longer quadratic in the canonical variables and one could thus 
expect strong quantum back-reaction even in the free model. However, this turns out not 
to be the case since one can choose non-canonical variables p together with J = pe*^ in 
which the system becomes linear: We have classical Poisson relations 

{P, Jjdsss = -iJ , {P, ^}class = iJ , {J, ^}class = "^ip ■ (21) 

Moreover, the Hamiltonian H = Yi^J — J) is a linear combination of the basic variables, 
making the dynamical system (p, J, H) a linear one. For such systems, if they remain 
linear after quantization, no quantum back-reaction occurs. (A quadratic Hamiltonian in 
canonical variables is a special case of a linear system.) That the system does remain linear 
after quantization has been demonstrated in |14j . 

In the presence of a potential the Hamiltonian become^ 

H = -P'n<P), (22) 

^As written, this refers to exponentials depending only on connection components while recent investi- 
gations of loop quantum cosmology have suggested the use of triad dependent holonomies, i.e. a{p) in e*"'^, 
in Hamiltonian constraint operators [121 [22l [23] . While basic holonomies in the holonomy-flux algebra 
do not depend on triad components, the appearance of triad dependent holonomies can be motivated by 
lattice refinements of an inhomogeneous state occurring in a physical state. Our following discussion of 
the qualitative behavior does not depend much on which form of holonomies is used in the constraint since 
alternative cases can be mapped into each other by canonical transformations of (p, J) before quantization. 
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which we again treat perturbatively. The expansion in V reads 

H = - ^j^V{ct>) + 0{V\<P)). (23) 

After quantization, we again obtain an equation of the form (1191) . In the variables used 
here, (p, J), we can use loop techniques and quantize J using holonomies. If we choose the 
ordering J = pe*^, basic commutators are 

[p, J] = hJ , [p, J^] = -h.P , [J, .P] = -2h{p + h/2) . (24) 

Quantum Poisson brackets obtained by taking expectation values of these equations are 
thus corrected compared to the classical ones: 

{p,J} = -iJ , {p,J} = iJ , {J, J} = 2t{p + h/2) . (25) 

Moreover, in these variables there are more than three fluctuation parameters, namely G^^, 
QJJ ^ QpJ ^ QP-^ , G^^ and G^^ , since we are using one complex variable J. These quantum 
variables are defined analogously as with canonical variables, e.g. 

G'^^"= l(Jjt + jtj) _ jj. 

They are restricted to the correct number by reality conditions G^^ = G^-^, G*^"^ = G-^-^ as 
well as 

-{P + W = - ^'^'^ + = - ci • (26) 

as a consequence of JJ'^ = p^. (G'^"^ — G^^ =: Ci turns out to be a constant of motion for 
the free model. Due to its appearance in the reality condition it plays an important role 
which will also be seen in the perturbed models discussed in what follows.) 

Dynamically, the difference to the Wheeler-DeWitt case is that H m a. metric or triad 
representation is a difference rather than differential operator [SUED]. Instead of a partial 
differential equation in the presence of a scalar we then have a partial difference-differential 
equation (see also [251 [26]). This is recognizable in our treatment only indirectly since we 
only implicitly deal with constraint equations for states. We rather express the information 
contained in the equation for a state through equations implied for expectation values, 
fluctuations and, if needed, higher moments of the state. The main effect then is the 
presence of J or e*'^, not c itself, in equations of motion. Also effective equations will then 
initially depend on e**^ instead of c, which suggests that they are simply obtained by a 
replacement of c by sin c as the real quantity having c as the low curvature limit for c ^ 1. 
However, this replacement on its own overlooks coupling terms of fluctuations which will 
be analyzed in this paper to provide a complete set of effective equations. 

The replacement of c by sine is complete for the free model as shown in [14J and 
indicated by the numerical studies of [Tni[II]. In this case, the quantum Hamiltonian lacks 
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coupling terms from fluctuations and free solutions 



Pv=o(0) = \{Ae-^ + Be^)-h (27) 
Jv=M) = \{Ae-^-Be^)+iH, (28) 

for expectation values and 

= \{c,e-'^ + c,e''')-\{ci + c,) (29) 

1 



G'vU{(l>) = -(c3e-2<^ + C4e2'^) + ^(3c2-Ci)-2z(c5e^-C6e-'^) = G^^io(0) (30) 
1 



G^yLo(0) = i^{c^e-^'^ - c,e^^) + ^{c,e^ + c,e-^) = Gl'Uct>) (31) 
G'Z,{ct>) = licse-'U c,e'*) + ^i3c^ - c) . (32) 
for all fluctuations can be determined explicitly. 

3 Effective equations 

Hamiltonian operators determine the evolution equation (|T9l) for a wave function ip. Equiv- 
alently, a state can be described in terms of its moments, or its expectation values of basic 
operators together with quantum variables. As already used in (|T0|) and seen in detail 
later, a Hamiltonian operator can then also be used to derive equations of motion for these 
infinitely many variables directly, without taking the detour of a wave function. These in- 
finitely many equations are in general coupled and difficult to analyze. Effective equations 
provide a well-defined scheme to control these equations by formulating them as equations 
of motion for the classical variables c and p identified with expectation values, or for a 
larger but finite set including also some of the quantum variables. Such equations can be 
obtained by systematic approximations (mainly the semiclassical one) to neglect most of 
the moments, or by solving for all but finitely many of the quantum variables in terms of 
those variables to be retained and inserting those solutions into the equations of motion 
for the variables to be retained. In the simplest case, one thus obtains effective equations 
for c and p amended by terms such as an effective potential depending on c and p. This 
is different from other methods such as perturbation schemes where solutions are already 
obtained as functions of time; one would solve all the equations up to a certain order but 
including those for c and p to find solutions. Solving for or eliminating only quantum 
variables but keeping a finite set of variables free requires different techniques. The result 
of such a procedure, a set of effective equations, will be much more powerful than pertur- 
bative solutions since it will illustrate quantum effects more directly and more generally 
by equations of motion (e.g. containing an effective potential) rather than only through 
specific solutions. 
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In general, the derivation of effective equations is more complicated than perturbative 
solution procedures, and it is not even guaranteed that effective equations of a given type 
exist to describe a certain regime of interest well. The key step in their derivation is a 
decoupling of the infinitely many equations for moments of a quantum system such that 
almost all of them can be solved at least approximately (or shown to be negligible) in terms 
of a finite set of remaining variables. The main tool to achieve this is often an adiabatic 
expansion for quantum variables which, for instance, provides the low energy effective 
action for anharmonic oscillators in quantum mechanics. But quantum variables such as 
fluctuations need not behave adiabatically in a given regime, or only certain combinations 
may behave so. Thus, in general the procedure of effective equations is not guaranteed 
to be of success. As we will see now, quantum cosmology provides examples for regimes 
where effective equations of adiabatic type do not exist. Nevertheless, as we will also show, 
bounces in loop quantum cosmology can, fortunately, be described by effective equations 
even though they turn out to be rather involved. 

3.1 Effective equations in quantum cosmology 

As shown in [T^ [15] and recalled in the preceding section, for a specific ordering of Hamil- 
tonian constraint operators the equations of motion ( ITOl) for expectation values of the triad 
p and connection c describing an isotropic universe sourced by a free, massless scalar do not 
depend on any quantum variables. They can not only be taken immediately as effective 
equations without genuine quantum corrections but even describe the system exactly. This 
is analogous to the equations governing expectation values of a harmonic oscillator or free 
field theory state. It is true not only in the Wheeler-DeWitt quantization but even in a 
loop quantization of the same model in appropriate variables. The replacement of connec- 
tion components c by "holonomies" sine in the constraint and the Friedmann equation, 
as often done as a shortcut to loop effective equations, is thus strictly justified for these 
models. In the presence of a potential, however, the situation is more complicated and ad- 
ditional quantum corrections arise from the back-reaction of fluctuations on the behavior of 
expectation values. These corrections can be crucial and must be determined for a robust 
analysis of quantum systems. This is analogous to adding an anharmonic contribution to 
the harmonic oscillator, or interaction terms to a free quantum field theory. 

In the latter, well-known cases an adiabatic approximation can be employed usefully to 
describe the quantum behavior around the ground or vacuum state of the system. Some of 
these corrections take into account properties of the interacting vacuum state as opposed 
to the free vacuum, or of dynamical coherent states as opposed to just kinematical ones. 
Quantum corrections can be computed order by order in an expansion in h combined with 
the adiabatic one. While the ^-expansion is always of interest to capture quantum effects 
in semiclassical regimes, an adiabatic expansion has to be justified independently. For 
perturbations around ground states of the harmonic oscillator or free field theories this is 
clearly reasonable: The exact states (simple Gaussians) of the unperturbed theories have 
constant fluctuations, and thus perturbations should result only in a weak time dependence. 
This expectation is borne out self-consistently by the fact that the adiabatic approximation 
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even with interaction terms provides well-defined solutions. 

The situation is different, however, for the quantum cosmological models studied here. 
We do have a solvable system which can provide the zeroth order for such an expansion, 
but it lacks a ground state or an otherwise distinguished state in which fluctuations and 
other quantum variables would be constant. Instead, exact solutions (|T2l) -(fT6l ) for the 
Wheeler-DeWitt model and fl27j) - fl32l) for the loop quantization in the absence of interac- 
tions show that most quantum variables are exponential functions of time. It is thus not 
clear a priori which combination of quantum variables, if any, can be treated adiabatically. 
Nevertheless, any candidate can be analyzed self-consistently by trying to set up an adia- 
batic approximation. As we will see in the descriptions that follow, several subtleties arise 
in this process. It clearly demonstrates that for general systems in loop quantum cosmol- 
ogy there is much more to an effective analysis than a simple replacement of connection 
components c by sin c. We will discuss this issue further in our Conclusions section, after 
providing a detailed analysis of effective descriptions. 

3.2 The Wheeler-DeWitt model 

We start with an attempt to develop effective equations based on an adiabaticity assump- 
tion for quantum variables. While this will ultimately be unsuccessful, it illustrates the 
main procedure and several difficulties which can arise. After that we describe possible 
choices of other combinations of the variables which may be assumed to be adiabatic, al- 
though none of them will result in well-defined effective equations for the Wheeler-DeWitt 
model. In view of a comparison with the loop model, and to confirm that the analyzed cases 
do not lead to consistent effective equations, it is nevertheless instructive to go through 
some of the details. In the next section, we will study a different, higher- dimensional type 
of effective equations which does exist. 

To derive effective equations from a Hamiltonian operator H one first writes the ex- 
pectation value of if in a general state as a functional of expectation values, spreads and 
higher moments, captured in quantum variables 

G'^'" = (((C-(C)) — (p-(p)r)weyl> (33) 

for 2 < n G Z and < a < n. This is possible because, with certain restrictions such 
as uncertainty relations discussed below, the quantum variables together with expecta- 
tion values allow one to reconstruct the state ■?/'(c, p, G"'") which has been used in com- 
puting the expectation value of H. Thus, the quantum Hamiltonian Hq{c,p,G"''"') = 
{iIj{c,p,G°'''^)\H\i/j{c,p,G"'''^)) is a function of expectation values c and p and quantum 
variables G"'". Due to the square root present in H it is difficult to write an explicit 
Hamiltonian operator in our case. Fortunately, this can be circumvented in our framework 
where we do not work with states and operators directly but with the coupled dynamics 
of expectation values and quantum variables. We thus implicitly define a Hamiltonian 
operator, order by order in a perturbation expansion as follows. 
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The quantum Hamiltonian can be obtained by formally expanding 
Hq = {H{c,p)) = {Hic+ic-c),p+ip-p))) 

n\\a) dp^dc''-" 

n=0 a=0 \ / -f" 

in c — c and p — p (with G""'^ understood as taking the value one, i.e. we use normalization 
of states; for n = 1 we have G"'^ = by definition). Higher order terms in these expan- 
sion parameters then combine to correction terms to the classical Hamiltonian containing 
quantum variables: 

Hq = cp- \^V{<P) - ^V(0)G^^ + (l + ^V^(0)) G^P - ^V{<P)GPP + ■■■ (35) 

where we restrict ourselves to positive p from now on. Since we treat (j) as an internal time 
variable, we do not include quantum variables of (f) and or any correlations between mat- 
ter and metric degrees of freedom. We do, however, include fluctuations and correlations 
of the gravitational variables relevant to the given orders, which is second in moments and 
first in the potential. Note that this suggests correction terms in an effective potential of 
the form 

Kff(0) = {l+c~^G'%<P)-2c~^p~^G'P{<P)+p~^GPPmyi.<P) ~ {l + C^-2C,+C,)V{cp) (36) 

(if we use the free solutions in the second step) and a "zero point energy" ~ C4, which 
is independent of the expectation values, to the Hamiltonian. To use this for an effective 
Hamiltonian in terms of c and p only, we would have to solve consistently for the quantum 
variables by appropriate means. A precise realization, as we will see, turns out to be more 
complicated. 

Equations of motion for expectation values and quantum variables are derived from 
quantum commutators determining Poisson brackets between expectation values and all 
the quantum variables. For instance, we have {c, p} = {[c, p]) / ih = 1 and, using the Leibniz 
rule, 

{G^^G^ = {{p')-{p)M{cp + pc)-{c){p)} = -2GPP (37) 
^QccQcpy = 2G^' (38) 
j^Qcc^QPP^ = 4G"P. (39) 

General formulas are provided in [16]. Equations of motion on this infinite dimensional 
phase space derived from the quantum Hamiltonian Hq{c,p,G"''"') = (H) agree with the 
quantum evolution equations determined from commutators with the Hamiltonian opera- 
tor. From the quantum Hamiltonian fl5S]) we then obtain equations of motion 

p 
c 



{c, Hq}=p-^Ji + c-'G'^^ - (cpr'G'^P) V(0) , 
p = {p,HQ} = -c-^^[l + ^c-'G^^-A{cp)-'G^P + y'GPP^V{<P), (40) 



Qa,2 ^ r^a,2^ ^ | ^ (0)G"^-i,2 + 2(1 - a) fl + U)] G""'^ - (2 - a) -V U)G''+^'^ 

\ / c 
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3.2.1 Non-adiabaticity of quantum variables 

To analyze possible adiabaticity of fluctuations, we use their equations of motion 



= 2 (l + - ^V(0)G^^ (41) 



^2 



QPv ^ 2^V{<t>)G'^ -2(l + ^V{<P))G^^ . (43) 

These will be the only equations needed for the order we are interested in. An adiabatic 
approximation proceeds by expanding the fluctuations as = Yl^=o A'^G"'^, replacing 
d/d0 by Ad/d0 in the equations of motion and expanding the resulting equations in A. 
One can then solve order by order for G^'^ and sum the contributions. The adiabatic 
approximation to a given order N is the sum G"^ = {J2e=o ^^^e'^)\>^=^- example 
of the low energy effective action for anharmonic oscillators shows, the zeroth order is 
already sufficient to derive an effective potential, and quantum corrections to the mass will 
be obtained at second order |T6l IT7] . 

Proceeding in this way, we obtain equations 

= 2 (l + ^\/(0)) G- - ^\/(0)GS^ (44) 

= ^Vi<P)Gl^--Vi<P)Gl\ (45) 

= 2^V(0)GS^ - 2 (l + ^r(0)) Gg^ (46) 

solved by 

= '''' 

Gff = ^Gf (48) 

in terms of Gq'^. Only these two relations follow since fim) - fHBl) presents a degenerate 
system of linear equations. This is a general property in these variables since the equations 
are derived with the degenerate Poisson algebra (137|) - (!39|) of three fluctuation variables. 
The free variable Gq'^, which in general can be a function on phase space and of 0, has to 
be determined at least up to a finite number of constants which specify the state used to 
expand around. Such a restriction arises by considering the next order: At first adiabatic 
order we have equations 

G^^ = 2{l + ^V{<P))Gr-^V{<P)GT, (49) 
Gl' = ^V-(0)Gr-ir(0)Gf , (50) 
Gl' = 2^\/(0)Gr - 2 (l + ^\/(0)) Gf (51) 
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which is again a degenerate hnear system, but now an inhomogeneous one for G^'^. This 
can only be solved for the G"'^ provided that the relation 

= ?!Gf + Gg. _ 2^ (l + GJ. 

is satisfied, a differential equation which in fact determines Gq'^ up to a constant. One can 
solve this equation for arbitrary potential V{(f)) by 

G- = , "^^'^^ (53) 

with the constant of integration x, and thus 

^cp ^ xc{l+pV{(t))/c^) 
° v/c2 + 2pV{<P) ' 

G- = ^£I(£V^. (55) 

Note that the initially worrisome inverse V{(f)) in (H7j) drops out of the final expressions 
which are perturbative in V{(f)). Nevertheless, having to divide by the potential in an 
intermediate step will turn out to be problematic. This feature will be cured in the loop 
quantized model. 

These solutions can now be inserted into the equations of motion for c and p to obtain 
candidate effective equations to zeroth adiabatic order 




X 



p^c^ + 2pV{(P) 



- + ^V{<j>))V{<P) + 0{{pV/c'f) (56) 



C3 

p = -p + -^\i -J^^=]V{<P) 

' Py^^+Wi^ 

1 _ 1^ + ^Vi<P)) V{<P) + OiipV/c^f) (57) 
cp J 

There is no h in these effective equations, but the free variable x parameterizes fluctua- 
tions of a state. It should thus be expected to be of the order h for effective equations 
corresponding to states saturating the uncertainty relation 
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However, when trying to implement these relations, we recognize that the effective equa- 
tions are inconsistent because there is in fact no state corresponding to this evolution: 
uncertainty relations are not satisfied. Inserting our solutions to zeroth order in the adia- 
batic expansion, we obtain 

G-Gl^ - {C^i? = -y^, (l + {Gir < 

for pV <^ . Thus, there is no choice for the only free parameter x in G'{^ which would 
make the adiabatic solution consistent with the uncertainty relation. This illustrates how 
one can self-consistently test whether an assumption of adiabaticity for certain variables 
results in reliable effective equations. Although the adiabaticity assumption for quantum 
variables has not resulted in consistent effective equations, Eqs. (156|) and (1571) suggest that 
true effective equations are more involved than just implying an effective potential (1361) : 
There are new expectation value dependent coefficients, and higher powers of the classical 
potential. In what follows we will look at other adiabaticity assumptions, before turning 
to higher dimensional effective equations in the next section. 



3.2.2 Non-perturbative potential 

Before changing our adiabaticity assumption, we probe whether the perturbative treatment 
of the potential may be responsible for the non-existence of effective equations. There 
could, after all, be regimes in which the potential is important in making quantum variables 
adiabatic. We thus start from a non-expanded classical Hamiltonian H = p^ (? — pV{(f)) 
and obtain the quantum Hamiltonian 

H- r2 — iTz:^ r^cc , - cpV{(f))/2 



2(c2 -py(0))3/2 

up to second order in quantum variables. Now, the zeroth order adiabatic solutions take 
the form 

_ c c^-pV{<j))/2 
^0 - V{^)c^-3pV{c^)/4' 

2 

= ^^wim^' ^^^^ 

and the uncertainty product 

is still negative independently of the potential in the allowed range \pV{(l))\/c^ < 1. (Allow- 
ing larger potentials would give only a narrow range 1 < pV{(f))/c^ < 2/^/3 ~ 1.155 where 
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the previous expression would be positive. But then no real is possible. Thus, drop- 
ping the assumption of a perturbative potential does not, by itself, give rise to consistent 
effective equations. 



3.2.3 Non-adiabaticity relative to expectation values 

While none of the quantum variables in (fT2l) -(fT6l ) is constant for a given state of the free 
Wheeler-DeWitt model which one may choose to perturb around, relative fluctuations 
G'^'^ !(? and G^^ /p^ are constant, and similarly for other quantum variables. One may thus 
try an adiabatic approximation for those combinations rather than one directly in terms 
of quantum variables which fails, as we saw. Let us now assume that variables of the form 

behave adiabatically where / is some fixed integer to be chosen. For the free solutions, the 
ga.,n ^jjj constant for any / since cp is constant, and thus adiabaticity seems more likely. 
(The constants Cj in the free solutions are not necessary to keep in the definition of g""'"^ 
since they can be absorbed without changing adiabaticity properties.) 

Equations of motion for the new variables follow from the product rule: 

_ PV{4>) /^2^cc I 3 1 2_ _ 2_o^T//j,\//iUPP 



where we ignore quadratic terms in gf"'". The zeroth order adiabatic solutions are 



CP _ 3 (/ + 8/3y - (/ + 2)p\/(0) 

^° " 4 c^-3pV{<P)/4 ^° ' ^^^^ 

PP ^ c\{l + 8/3y- {I + 2)pV{4>)) 

^° (c2-3py(0)/4)((/-2/3)c2-/p\/(0))^o ^ ^ 

which solves the remaining zeroth order equation only for / = 0, / = 1 or / = 2. 

The case / = 0, i.e. an assumption of adiabaticity for (7^^ = G^^/c^, g^^ = G'^^/cp and 
gPP — Qvv jjP' ^ gives solutions c/g^ = 2g'^ and g^ = Ag^^ which in this case are actually exact 
to any adiabatic order. However, the uncertainty product c^p'^g'^'^g^^—{cpg^'''Y now vanishes 
for any choice of free parameters and uncertainty relations are still violated. This happens 
only marginally since for a slightly positive value one could choose the free parameter, 
analogous to x before, large enough to ensure that the uncertainty relation is satisfied. It 
could happen that higher orders of the adiabatic expansion would make the uncertainty 
product positive, but we do not pursue this here since, if successful, it would most likely 
give only special regimes where effective equations would be valid. 
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For / = 1 we obtain 



with a negative uncertainty product. For / = 2, we have 

with a vanishing uncertainty product. None of these cases gives a consistent set of effective 
equations. 



3.2.4 Additive non-adiabaticity 

Instead of factorizing off the time dependence of the free solutions as in the variables 
ga,n^ one can try to subtract it off, i.e. attempt an adiabatic expansion in G'^'^ — C^c^, 
Qcp _ g^j-^fj Qvp _ C^p^. In contrast to the previous case the equations now also 

depend on parameters Cj of the free solutions. It is easy to derive the resulting equations 
which, at any adiabatic order, provides a set of inhomogeneous linear equations for the 
quantum variables. There are three equations for three fluctuation parameters, but with a 
vanishing determinant of the linear system. It turns out that the inhomogeneity resulting 
from subtracting off the zeroth order solutions does not allow any solution. Thus, also 
adiabaticity additive to the zeroth order solutions is not realized. 



3.2.5 Non-adiabaticity with respect to free solutions 

Had one of the previous attempts succeeded in providing effective equations, the result 
would have been time dependent only through the scalar potential V{(j)) as is clear from 
the examples provided. Our final attempt brings in even stronger time dependence by 
assuming again adiabaticity relative to free solutions, but using the free solutions (fT2l) " 
( IT6l) in explicit, time dependent form Gy'^Q^ct)) rather than implicitly as, e.g., Gy=o 
before. Thus, we now assume that variables := G''^'"'/G'y"Q(0) are adiabatic. The new 
variables satisfy equations of motion 



■a,n ^ _ G^y=o(0) a,n 



or explicitly 



^?^,^.o(0)^cc _ V^^)^^MM±^My. , (66) 



fc2 
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To zeroth order, adiabatic solutions would be given by 



which, when inserted into equations of motion for c and p would give time dependent terms 
not only through V{(j)) but also through Gy'^Q{(f)). This is an exponential time dependence 
and thus stronger than that of a polynomial potential V{(j)). However, also here the result 
is not consistent since the uncertainty product 

to leading order in pV{(f))/c'^ is again negative. 

We did not manage to find a choice of adiabatic variables which would give rise to 
consistent effective equations. However, since some cases led only to marginal violations of 
the uncertainty relation, one may still be hopeful that a choice exists. More importantly, 
changes to the equations of motion such as those by a loop quantization could possibly 
lead to the existence of consistent effective equations which we will probe next. 



3.3 Effective equations of loop quantum cosmology 

The solvable loop model is formulated in non-canonical variables, satisfying a Poisson 
algebra different from that of the Wheeler-DeWitt model. Thus, also the procedure of 
deriving effective equations is different and has to be analyzed anew. We must now consider 
two different sets of quantum variables, each one related to J or J. We will first obtain 
general equations of motion for these variables, which are then to be restricted by reality 
conditions ensuring that J J = p^ or, at the quantum level, J.P = fP'. Expanded in powers 
of quantum variables, the quantum Hamiltonian reads 

^ 2i J-J {J-J)^\ J 

pMfL _ ^p A _ 3,PmG- + . . . . (68) 

From now on we ignore terms of higher moments since they will be subdominant in a 
semiclassical expansion. Higher moments of typical semiclassical states are smaller at 
higher orders. Moreover, a moment of order n appears in the quantum Hamiltonian with 
a coefficient obtained by an ra-th derivative of the classical Hamiltonian; see flM|) . Higher 
moments of order n are thus suppressed by additional small factors of the form p~^{J — 
j^-n+k compared to the classical perturbation term. There is a subtlety with this argument 
when looking at equations of motion: In Poisson brackets {G'^'", G'*'"^} with n + m > 4 
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(formally including n = 1 or m = 1 as expectation values) there can not only be linear terms 
of moments of order n+m—1, as in the basic relations (!69!) - (!771) below, but also terms of the 
form (^.^'(^c,n+m-i-« 0<Z<'n, + m — 3, as illustrated by the examples in App. |Al (Also 
h'^-terms appear, but they play a role only at higher orders.) Depending on the powers I 
realized, the expectation value dependent prefactors could, in equations of motion, cancel 
suppressions involved in coefficients of higher moment terms in the Hamiltonian. This is 
potentially problematic if the highest allowed powers for / do occur. For instance, we have 
{G"'''^, G^'^} of the form (■)G"^'^ with a single factor of expectation values, but in a bracket 
like there could be terms {■)^G'^'^ with up to two factors. Then, the additional 

suppression of G*'^ in the Hamiltonian would be cancelled in equations of motion by the 
additional factor in Poisson brackets. It would require one to take into account all terms 
containing C^'" in the quantum Hamiltonian, which would clearly be unmanageable. 

Fortunately, this problem does not arise and there is a consistent way to derive a trun- 
cation for equations of motion to any given order: For equations of motion for expectation 
values, Poisson brackets only have terms linear in moments G'''"' of the same 

order as used in the bracket, and possibly terms {■)G'^'^~^ with a single factor of an ex- 
pectation value. (There may also be terms with higher powers of h which we can safely 
ignore.) The terms with a single factor of expectation values may arise because taking a 
commutator of, say, p with (J — J)" as it is needed for {p, G"^"} replaces one factor of J — J 
by just J itself. To bring the result into the form of a moment of order n, one has to add 
J multiplied with a moment of order n — 1. The only other terms that may arise in such 
cases are due to reordering, which contributes higher powers of h. (They would be relevant 
at higher moment orders than considered here.) Taking into account suppression factors 
in the Hamiltonian, equations of motion for expectation values thus have suppressions of 
at least (•)^~"' in higher moment contributions. Being n-dependent, this leads to higher 
suppression for higher moments and provides a consistent truncation. The same order of 
truncation can then be extended to quantum variables by requiring that the reality con- 
dition fl26|) is preserved by the truncated equations of motion. This consistent truncation 
procedure is illustrated below for moments of up to second order. 

While all quantum variables commute with expectation values when defined in terms 
of canonical variables, quantum variables and expectation values do not Poisson commute 
for non-canonical variables such as p and J. Following the general scheme used before, 
Poisson bracket relations between the variables are derived as 

{G'',p} = = {GPP,p} (69) 
{G'f^,p} = iGP-^ , {GP-^,p} = -iGP^~ (70) 
{G-^-^,p} = 2iG-^^ , {G-^~-^~,p} = -2iG-^^~ (71) 
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and 



{GPP, J} = -2iGP^ , {GPP, J} = 2iGP^ , (72) 

{G^-^,J} = -2iGP^ , {G^^,J} = 2iGP\ (73) 

{GP-^,J} = -iG^^ , {GP-^,J} = iG^^ + 2iGPP, (74) 

{GP^,J} = -iG^^~-2iGPP , {GP^,J} = iG^~\ (75) 

{G^-^, J} = , {G-^-^, J} = AiGP^ , (76) 

{G^'-^, J} = -AiGP^~ , {G^^, J} = (77) 



between quantum variables and expectation values. Notice that fl69l) . fl72|) and fl73l) imply 
that the difference G'^^ — Gpp commutes with any function of p, J and J only. 

Also in contrast to the Wheeler-DeWitt model whose canonical variables gave rise to 
a closed Poisson algebra for the fluctuations C^'^, Poisson brackets between fluctuations in 
the loop model involve higher moments due to the non-canonical form of the basic variables 
{p, J) . Poisson brackets between two fluctuations thus in general depend on a quantum 
variable of order three, but they also involve fluctuation dependent terms. While the third 
order variables can be cut off in a semiclassical treatment to the order we are interested in, 
the latter terms, as we will see below, have to be included in a consistent set of effective 
equations. We thus need to determine the full set of Poisson brackets between quantum 



variables of order two. As derived in Appendix \^ they are 

{GPP,G^^} = -2i{JGP^~ - JGP-') (78) 

{GPP,G-^-^} = -AiGP-^^ - AiJGP-^ = {GPP, G^^} (79) 

{GPP,GP-'} = -2iGPP-^ - 2iJGPP + \ih^J = {GPP, Gf^"} (80) 

{G^-\G'^-^} = -MGP^^ - Ai{p + h/2)G-^'^ = {G^^, G^'^'} (81) 

{G^'^,GP^} = -2iGPP^-2i{p + h/2)GP^ +iJG^^~ -iJG^^ +\ih^J = {G^^, Gp^} (82) 

{G-^^,GP-'} = 2iG^^-^ + 2iJG^-^ = {G-^J, Gp^} (83) 

{G^-\GP^} = GiGPP'^ + 8i{p + h/2)GP^ -2iJG-^-^ -ih^J ={GJJ~G^ (84) 

{G'^^,G'^^} = SiGPPP + Wiip + h/2)GPP + 8iip + h/2)G^'^ ~8iJGP^~ -SiJGP'^ (85) 
+4iffp + 2ih^ 

{GP^,GP-^} = AiGPPP + Qi{p + h/2)GPP -iJGP^~ -iJGP^ + \ih^p+\ih\ (86) 



The relevant relations for leading order back-reaction effects in the loop model are brackets 
of the form {C^'^, (■)} together with contributions of {G**'^, G'''^} from moments up to 
second order. (The fluctuation independent ^^-terms can to this order be ignored as well 
since G"'^ are typically of order h near saturation of the uncertainty relation.) Expectation 
values and second moments thus provide a closed system under taking Poisson brackets in 
this form, albeit a non-linear one. 
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3.3.1 Equations of motion 

Just as the algebra ( !37l) - (!39l) of fluctuations in canonical variables is degenerate, being 
given by a Poisson structure on a 3-dimensional manifold, the relations {G"'^, (■)} are 
degenerate in the sense that functions of the fluctuations exist which commute with both 
expectation values p and J. It is easy to see that, as already noted, one such function is 
linear, G"^'^ — = Ci. The second function is quadratic, 

C ■= 2(G^'P)2 + {G-^-^f - + . (87) 

(Its value is G = |c^ — iciC2 + |c2+4c3C4+4c5C6 for the free solutions fl271) - fl32l) .) The linear 
function has one immediate implication: It is exactly this combination which appears in 
the reality condition ( 126|) implied by J J'^ = in the loop quantization. In the free model, 
the combination ci = G'^^ — G^^ is a constant of motion. Whether or not this remains true 
when perturbations are included is an interesting question because the form of the reality 
condition (corresponding to an implementation of the correct physical inner product) is 
essential for singularity resolution: it picks out exactly the bouncing solutions in the free 
model (unless the constant of motion ci would be negative and |ci| large, which is never 
realized for states of the free model being semiclassical at least once). If the form of the 
reality condition could be shown to be preserved by quantum corrections, there awould be 
no obvious violations of the bounce as one could in general imagine if ci would not remain 
constant and could evolve to large negative values from initial semiclassical states |15]. If, 
however, G'^'^ — G^^ does not remain constant when the free model is perturbed, one would 
have to be more careful with conclusions about the bounce. 

Since the difference G'^'^ — G^^ commutes with any function of the expectation values, 
its preservation does not depend on the precise form of the classical Hamiltonian. However, 
for leading quantum corrections we have to consider fluctuation dependent terms in the 
quantum Hamiltonian ( !68l) . which contribute extra terms to the equations of motion of 
fluctuations. These fluctuation dependent terms do not necessarily commute with G*^"^ — 
GPP. 

Using the Poisson relations, we obtain 

P = -^ + j~-,p'V{<P) + 3jj±d^y{G''^ (88) 
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and 



j ^ -(p^fi/2)+(j^ + 2f^)m (89) 

2 

-3— ^^(G-^-^ - G"^"^' - 2G'^'^')y(0) + 6-^Gf"^y(0) 
(J — J)^ J — J 

for the expectation values and 

G^" = - (l - (J? 7)2 ) + G^^>('/') 



4p 



.3 



[jG^''^ - JGP-^ - JG^^ + JG^^) 1^(0) 



(J - J) 



,6.^(^+^Vy(0)-^^i^n0) (90) 



(J-J)2 2 {J-Jf 



and 



Q„2 _ _ __ 

-,J (J + JlG-^"^ - (2p + mfG^*^ + G^-^) - JG^^ - JG^^ 
{J -J) 

+\h\J + J))V{<t^) + ~{JG^' - JG^')V{ci>) 
^ -(l-3p- ^^^(^yiy,^/^)) )(G^- + G^--)y(0) 

_3P^^^,. _ ^..^^^^^ - (1 + iG^^ - 
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for the fluctuations featuring in the reahty condition, ignoring terms quadratic in fluctua- 
tions. A combination gives 



QjJ _ Qpp 




(92) 



which can be seen to equal 2d(| — {p + /i/2)^)/d0 from (!88|) and (!89|) . Terms included 
in these equations of motion preserve the reality condition. Thus, this consistent trunca- 
tion of the quantum equations describes the physical evolution of expectation values and 
fluctuations as they would be computed from a physical state. 

However, these equations also show that the contribution G^^ — G^^ to the reality con- 
dition is no longer constant when fluctuation dependent corrections in Hq are allowed for. 
It is then possible that the term ci, which was constant in the free model and determines 
whether or not expectation values bounce, changes in time so as to prevent a bounce when 
small volume is approached during long evolution. A firm conclusion about the presence of 
a bounce for a massive or interacting scalar thus requires a more detailed analysis taking 
into account the behavior of fluctuations. We will come back to this issue in the next 
section. 

Equations of motion for the remaining second moments are 



^^^-^{-2{J + J)G^^ - 2JG^'' + 2{p + h/2){AG^'' - 2G^^) + 2JG^ 
[J -Jy 

h'j + \h\j) + ^y^^(6(p + h/2)G^^ - JG^' - JGP' + \h^p + 1^3) 




JJ 



(93) 



GpJ 
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and 

G'' = -2 (l - ^j^) G^'m + j^G^M^) (94) 



(J - J)3 



[8{p + h/2){2GPP + G^^ - G^^) - 8JGP^ - SJG^-^ + Affp + 2H^) 



[J — jy J — J 



= GJJ . 

As in the Wheeler-DeWitt context, we will now briefly describe attempts to find adiabatic 
regimes of quantum variables, followed in the next section by a perturbative analysis of 
the equations of motion. 

3.3.2 Non-adiabaticity of quantum variables 

For simplicity, we now take into account the classical Hamiltonian as a function of expec- 
tation values only in order to compute the leading order terms in equations of motion for 
fluctuations. These equations are 

G^P = -(l- y%,, V{(l))] {G''^ + G^^') = G^^ , (95) 



V {J -J? / 

G^' = -1 fl-— ^^V(0)Vg-'^ + G^-^" + 2Gpp) + ^^V(0)GP'^ = 6^,(96) 
2 \ [J — J) ) J — J 

G'' = -2(^l-^^^V{cl^)^G^' + j^V{cl^)G''=WJ. (97) 

These are not the complete equations of motion; we use them only to illustrate obstacles 
to the existence of adiabatic effective equations. The full equations of motion to this order 
will be used later on for a new type of effective system. 

To zeroth order in an adiabatic expansion of the fluctuations, Eq. (p5|) shows that Gq"^ 
is purely imaginary, from which Eq. (!97|) then shows that Gq'^ is real. Moreover, Eq. (!97|) 
gives the relation between these two quantities, 

° ° " l-2p3V(0)/(J- J)2 J -J ^^^^ 

which when used in Eq. (196!) leaves only one free function such as G^: to leading order in 
p^V {(p) / {J — J) we have 

Gi' + Gf + 6G- + 2ci = ^-^j^iGi' + ) (99) 

and thus 

^° -~l-18p4V(0)V(J-J)2- ^ ^ 
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In these equations, also Gq — = c\ has been used. Unhke in the Wheeler-DeWitt 
case, the remaining free function Gq^ would not be fixed at the next order in the adiabatic 
expansion but through the second, quadratic constant C in ( |87ll . 

We first verify that the fluctuations are at least positive as they should by definition. 
While this does not apply to G^'\ J not being self-adjoint, the quantity QJ+J'J+-J ■= 
((J + J'^Y) — {J + j"'')^ must satisfy positivity conditions. Eq. fllOOp in fact shows that 
Gq"'+Go'^ is negative, and thus Go'^+Gq'^+6Go^+2ci is positive from fl99l) . This is relevant 
because the latter quantity can be written as G"^^ + G^^ + QG^^ + 2ci = G-^+*^'"^+*^ + AG^p 
in terms of the positive fluctuation Q'^^'^''^^'^ := {{J + J'^Y) — {J + J^)^ of the self-adjoint 
J + The zeroth order adiabatic solutions are thus not in conflict with positivity of 
fluctuations. 

Next, we check uncertainty relations. This time, there are several independent ones [15] 
because we are using partially complex basic variables. Independent uncertainty relations 
must then be fulfilled for each pair (p, J + J"'"), {p,i{J — J^) and J + ,i{J — J^) of 
self-adjoint operators. These relations, derived in [TJ], take the form 

QPPQj+j,j+j _ ^QP,j+JY > -1^2(7-7)2 (101) 

QPPQi{J-J),^{J-J) _ (^qpMJ-J)^ > ^fi^^J^Jf (102) 

_ ^^J+J,^iJ-J-)y > h\2p + hf. (103) 

The first one is easy to satisfy to this order since G^'^^^ = G^^ + G^^ vanishes. The left 
hand side of (110 II) is then positive as just shown, and one can find appropriate scalings (i.e. 
sufficiently large C) to ensure the uncertainty relation noting that J — J is nearly constant 
for perturbative solutions. The second relation is more difficult to verify since G^'*^"^""^^ is 
non-zero. For the zeroth order solutions, we have 

^pp^i{j-J),i{j-j) I /-ip,i{J-J)\2 _ r'PPinJJ _L nJJ onJJ\ _l inv-^ r'pJ\2 

which is at least positive since we saw that Gq*^ + G^"^ is negative and (Gq"^ — Gq'^)^, 
although negative, is of order l^((/>)^ as a consequence of fl98|) . Moreover, using the free 
solutions as a guide indicates that one might be able to satisfy the uncertainty relation: In 
the free case, we have 

-(G^io + Gilo - 2G^io) = 2(ci - c^) = A{AHy 

which is positive and constant. With the factor of G^^, which in the free case behaves as 
at least far away from the bounce one may have a chance to satisfy f ll02p . Similarly, 
the third relation, which is simpler due to q-^+-^''^(-^~-^^ = "^{Gq'^ — Gq'^) = 0, has a left hand 
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side 



,J+J,J+J ^i{J-J),i' 

"-^o 



(J- J) 



J+J,i{J-J)\2 



{Gi' + Gi'Y + A{Gi'f 
(6G^^ + 2ci)2 



(1 - 18p4\/(0)V(J - J)2) 



-+4(Gr + ci)^ 



32(G^P)2 - 16ciG^P - 144 




for the first two leading orders in p'^V{(j))/{J — J)- It can be positive for states satisfying 
Ci < — 2Gq^, but also has to be larger than the 0-dependent Ah^p^ as required by the 
uncertainty relation. A look at the free solutions also here indicates that this might be 
possible: both ci and GyLo ^^^^ order h for saturation of the free uncertainty relations, 
and G^^Q behaves as far away from the bounce. 

Using the candidate solutions for adiabatic quantum variables in the equation of motion 
for p then suggests an effective equation 



up to higher order terms in the potential. Here, Ci is a constant parameterizing the state 
around which the effective equation is derived. Similarly, Gq^ depends on ci and the second 
constant G in (IHTll which is preserved by the equations of motion we consider here. Thus, 
we would have an effective equation containing only two state dependent parameters but no 
free functions. However, in the final step of determining Gg^(ci, G) an inconsistency arises: 
using the candidate adiabatic solutions in flHTj) for a constant G, G^ has to be constant up 
to perturbative terms of the order p^V"^ /{J — J)^. However, we saw that the uncertainty 
relations require G^ to behave as p^ as it does in the free case. With a potential, there is 
no adiabatic regime in which this can be satisfied and the effective equation fll04p is, after 
all, inconsistent. 

3.3.3 Non-adiabaticity relative to free solutions 

For the loop quantized model it is more complicated to introduce adiabaticity relative 
to free solutions because free quantum variables (|27l) - (|5^ are not simply proportional 
to powers of the expectation values as in the Wheeler-DeWitt case. The most direct 
possibility is to introduce the free solutions as explicitly time dependent functions Gy"Q(0). 
Resulting effective equations will then be strongly time dependent not just through the 
scalar potential. Then, for variables (7"'" := G"'"'/Gy"Q(0) and using (|95|) - (|971) we have 



p = 



-12 




) 



Gr\/(0) (104) 
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equations of motion 



^ \ / G^LM J Gfco(0) .A , G^yU<P) + Gtoi^) 

^ V {J - Jy J [G^yU<l>)^ G'yU<pf J G^yU<P) ^ ' 

^' (^-^)VlvG'^/=o(</>) G'^yio(</') G^iM 

g^io + g^io(0) + 2G^"=o(0) 6p^V(0) \ 

<=o(0) J-J 1^ ^ ' 




g'' = - 



J-J J 

( _ \ / Gfco(0) , gfco(0) „A , GfcoW + GfcoW ,j 

V iJ - jy ) \Gi/^,J ^GiiM ) GiiM ' • 

There is now a crucial difference to the situation in a Wheeler-DeWitt quantization as in 
Eqs. f l65p . (jnSD and fl^7|) : \^ = is not a singular point for the adiabatic approximation. In 
the Wheeler-DeWitt equations, some terms completely vanish in the free case such that the 
solutions change dramatically. For instance, solving for g^^ requires one to divide by V{(f)) 
which is not possible in the free case. Thus, the existence of consistent solutions satisfying 
uncertainty relations in the free case, and obviously being adiabatic for the variables = 
G"''"' / Gy^Q{(j)) = 1, does not guarantee the existence of consistent adiabatic solutions even 
under slight perturbations including a small potential. The equations obtained here for 
the loop quantization, on the other hand, are perfectly perturbative in V{(j)) and the free 
solutions do not present a singular point for the adiabatic expansion. 

Nevertheless, there are again no adiabatic solutions in the presence of a potential. 
While the linear system of six equations for six (real) variables at zeroth adiabatic order is 
degenerate when = 0, allowing a non-trivial solution where all g"'''^ = 1, the determinant 
of the coefficients can be seen to be non-vanishing for V ^ 0. Thus, the only solution in 
this case is g"''"' = which again violates the uncertainty principles. 

There might be adiabatic regimes for the equations (155]) - (PT]) and to this 

order, but this would be more involved and does not seem promising. We therefore turn 
to a new type of effective equations which does not require an adiabatic expansion. 



4 Higher dimensional effective equations and pertur- 
bative solutions 

Although our attempts to find adiabatic regimes and to solve for quantum variables in 
terms of expectation values were largely unsuccessful, prohibiting the derivation of effec- 
tive equations solely in terms of expectation values, effective equations for an enlarged 
system do exist. Rather than solving for the fluctuations, we keep them as independent 
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variables in addition to the expectation values. This is still a system of classical type, given 
in terms of finitely many variables, since we already removed higher moments as part of 
the semiclassical expansion. But it is certainly different from the exact classical system 
not only in the form of its equations of motion but also by the number of independent 
variables. From a general viewpoint, this is in agreement with expectations from effective 
actions which generically give rise to higher derivative terms. Higher derivative actions 
imply higher derivative equations of motion which, when taken at face value, describe a 
larger set of independent degrees of freedom: additional parameters have to be specified 
for a complete initial value problem of higher derivative equations. However, when higher 
derivative terms arise from a perturbation expansion not all the solutions are consistent 
within this scheme. Only solutions analytic in the perturbation parameter are consistent, 
which can be shown to allow a number of independent solutions agreeing with the unper- 
turbed case [271 [28]. Thus, higher derivative effective actions do not truly introduce new 
degrees of freedom. Keeping some of the quantum variables, on the other hand, requires 
additional initial conditions even for perturbative solutions; in this case new quantum 
degrees of freedom truly arise. 

4.1 Wheeler— DeWitt model 

In variables c and p, we thus have a set of five coupled effective equations 



Back-reaction terms of fluctuations are included in the equations of motion for expectation 
values, although the role of an effective potential is illustrated more indirectly than it would 
be for adiabatic effective equations. 

For perturbations around free solutions, assuming small values of the potential ^(0) = 
eV(0) with a small parameter e, we define 



and similarly for the other variables, where po is the zeroth order solution f|T3l) correspond- 
ing to V{(f)) = 0. (For simplicity we now use a simple zero as subscript for the free solutions 
rather than = 0" . Confusion will not arise because the adiabatic approximation will no 




(105) 



P(0) =Po(0) + epi(0) + ■ •• 



(106) 
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longer be used.) All free solutions are explicitly given by (fT2!) - (fT6l) . To first perturbative 
order in e, the equations of motion (11051) become 



= + I -- - 73 

Co Cq 



Go + \g^') V(0) 



PI = -Pi+(4§-S^o^ + 2|g--^G-)V(0) 



= 2Gf + 2 ( ^G^j" - -G^n V(0) 
V Cq Co j 

Gf = (^Gl' --G'i]v{<i)) 

V Co Co J 

V ^0 Cq J 

which are inhomogeneous differential equations whose homogeneous terms are linear. Pro- 
ceeding in the same way to higher orders, this provides a well-defined and self-consistent 
solution scheme for the interacting dynamics of the quantum system. There are two expan- 
sions of the full quantum equations: a semiclassical one (truncating higher moments) and 
the perturbative one in e (removing higher powers of the potential). When extending to 
higher orders, one would thus have to find out what the relative magnitudes of correction 
terms in the different expansions are. 

Using the free solutions ffT2]) - f[TB]) for the different variables we obtain the perturbative 
solutions up to first order in e: 



ci(0) = 5ie^+(C2/Ci + C2C3/Cf-C4M2)e^y V(r)e-='Mr 

pi(0) = Sse-^ - {-Cl/2Gl - 'iClC^l2C\ + ^C^C^ICl - ^5/2^2)6-'^ V(r)e-3Mr 
6-^(0) = B^e^-^ - 2{C2Cs/Cf - C,/Ci)e^'^ j'^ V(r)e-3-dr 
Gr{<P) = B,-{GlC,/Cl-C,/C,) j%{T)e-^^dr 
Gf(0) = B^e-^'^ -2{GlCi/Cl-C2C^/Cl)e-^^ j%{T)e-^^Ar 

where (5i,...,i?5) are further constants of integration simply adding up to the zeroth 
order constants in (fT2l) - (fT6l) . Specifically, for a mass term potential V(0) = 
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we have 



Pi{ 

GTi 
G'f( 



Bie^ + ^{C^/C, + C2C3/CI - C,/Cl){2 + 60 + 9(1)^-^^ 

526-* - ^i-Cl/2Cf - 3ClC3/2C^ + 2C2Ci/Cf - C^l2Cl){2 + 60 + 902)6" 

i?3e'* - ^{C^C-^ICl - C4/Ci)(2 + 60 + 902)6-^ 

- ^{ClC,/Cl - C,/C,){2 + 60 + 902)e-3* 
B^e-^'i' - ^{ClC/Cl - C2C,/Cl){2 + 60 + 902)e-5<^ . 



These solutions aheady display interesting properties of the system. Compared to the 
free solutions which are pure exponentials in 0, there are different types of corrections. 
Since we are perturbatively adding a potential one certainly expects corrections to the 
solutions which correspond to some of the polynomial terms. For the classical variables c 
and j9, these agree with corrections to the classical solutions resulting from the potential. 
But there are additional corrections displaying genuine quantum behavior since fluctuations 
Qa.,'2 evolve, too, and couple to the expectation values. Although they are structurally 
similar, one can easily disentangle both types of corrections by setting the zeroth order 
fluctuations to zero. (This is not possible for the quantum system since it would violate 
uncertainty relations. But we are free to do so formally to discuss different correction 
terms.) Zeroth order fluctuations are zero if we set C3 = C4 = C5 = 0. This removes some 
corrections from Ci and pi which can therefore be attributed unambiguously to quantum 
back-reaction. They differ from classical corrections only in their coefficients, related to 
the integration constants. The time dependence of the correction, in this model, is thus 
the same as the classical one and quantum back-reaction does not change the expectation 
values significantly. Moreover, their coefficients are small if a state is semiclassical at large 
volume, which requires that the fluctuation parameters C3, C4 and C5 are small compared 
to Ci and C2. 



4.2 Loop quantization 

The loop quantization, being formulated in non-canonical variables (p, J), leads to a dif- 
ferent algebra of expectation values and quantum variables. Moreover, due to the use of 
the complex variable J, which is typical for loop quantization based on holonomies, there 
are more independent variables and equations of motion. This is supplemented by a reality 
condition fl2B]) which leaves the correct number of degrees of freedom when implemented. 
The Poisson brackets fl7n]) - fl77|) together with the quantum Hamiltonian fIBSl) result in 
equations of motion ([HHD and ([89]). With equations of motion ([90]), ([H]), ([93]) and ([94]) for 
C^'^, this provides our higher dimensional effective system of the loop quantization. 

For perturbative solutions of the expectation values we expand as in (11061) . use free 
solutions Po(0); -^0(0) and ^0^(0) from (IT^ - f[T^ on the right hand side and integrate. We 
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do not require solutions or even equations of motion for G"'^ to this order in a perturbative 
solution scheme if we are only interested in pi{4>) to analyze the bounce. Strictly speaking, 
Gi^ and Gf'^ are also necessary to implement the reality condition, but we can avoid explicit 
solutions since we already verified that the reality condition is preserved under evolution 
generated by our truncated quantum Hamiltonian. The evolution of G{'^ — G^^, however, 
will be considered below to analyze effects of the reality condition on the bounce. 
For the real part of J which couples to p, the equations give 

liJ+jy = -(p + V2) + ^pV(0) + 2^^^^±i^V(0) (108) 



Q „2 __ _ 



2(J- J)2 

Zeroth order solutions fl27|) - fl32|) for this system now describe bouncing cosmologies. 
They satisfy Jq — Jq = 2iHQ which is the constant free Hamiltonian, and Gq-^+Gq"^— 2Gq'^ = 
2(c2 — Ci) (with Ci — C2 = 2(Aifo)^ > 0) which can be used to simplify some calculations. 
Perturbed equations of motion are more complicated than for the Wheeler-DeWitt model, 
but we can focus on a few dominant terms. Moreover, we can assume A = B in (1271) . 
( |28|) since we would just need to shift our internal time otherwise. We keep the classical 
perturbation, given by p^( J- J)"^y(0) for p and |p^\^((/)) + 2p^(p + ^/2)( J- J)"^V"(</)) for 
|(J+ J)', which contributes (considering only the highest positive and negative exponents) 



Jo + Jo 3 .^e-^^-e 
-Po A 



3 _ 

AHi ^ 32HI 

to pi and 

to |(Ji + Ji)'. This is to be compared with the quantum correction term having the 
strongest time dependence, which turns out to be 

for p and 
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for + J)' ■ The inhomogeneous linear differential equations for the first order pertur- 
bations are solved by 

= ^(A(0)e-* + 5i(0)e<^) 
+ Ji(0)) = i(Ai(0)e-*-i?i(0)e*) 



2' 2 



with 



for the perturbative potential again defined as V{(f)) = eV(0). Thus, 

+ \§2^i (e^"'^ f V{y)e^ydy - 6"^+^^ V{y)e-^ydy^ (109) 



1 ^ 

4^, 



with 



tti = a/c^ , "2 = \/c^ , e^^^ = — , e^^' = — . (110) 

C4 C5 

(Integration constants from integrating Ai{(f)) and -Bi(0) can be absorbed in the zeroth 
order ones.) The parameters 61 are related to the squeezing of the unperturbed state, with 
(5/ = for an unsqueezed state pTS] . 

To this perturbative order, the bounce persists even under inclusion of a matter po- 
tential. As in the free model, the key to deciding whether or not there is a bounce, i.e. 
whether or not p is bounded away from zero, is the reality condition (!26|) . In the free 
case, its implementation resulted in selecting the cosh-solution for p and ruling out the 
sinh solution, proving the bounce under the assumption of a semiclassical state at large 
volume. (The integration constants A and B must have equal sign in (1271) if Ci is not too 
large and negative.) Perturbatively, the reality condition becomes 

1/^2 ^ \J\2_(^p^^/2f_QPP^QJJ^^J^^2_^p^^^^2f^c^ (111) 

+2e (ReJoReJi + ImJoImJi - {po + h/2)pi - ^Gf + ^G^^") + 0[e^) . 

The zeroth order is already implemented by starting with the correct zeroth order solutions. 
The linear order can then be seen to be automatically satisfied as well, such that the 
bouncing behavior at zeroth order is not affected by the perturbation. 
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It is, however, also clear that the bounce, with its mixture of collapsing and expanding 
functions, leads to stronger quantum back-reaction than classically. Although ai and 0:2 
are small compared to for a state which starts out semiclassically, depending on the 
values of 5i and 82 there can be exponential magnifications. For C3 ~ C4 and C5 ~ cg, which 
corresponds to nearly unsqueezed states and implies fluctuations nearly symmetric around 
the bounce point |15] quantum corrections are still small at the bounce if they were small 
for an initial solution. But for squeezed states with C3 differing from C4 corrections can be 
noticeable even earlier, before the bounce. This is illustrated by the solution 

Gi'-Gl^ = 44f(ci-C2)((2 + 902)cosh(30)-60sinh(30)) (112) 
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1 A 

~108^ 



[c^ - ce) ((2 + 90^) sinh(30) - 60cosh(30)) 



- 1^ ai ((2 + 902) cosh(30 - 5i) - 60sinh(30 - 5i)) 

1 

+ — — as ((2 + 902) cosh(30 - ^2) - 60sinh(30 - ^2)) 
12 -no 

for the combination of fluctuations featuring in the reality condition (specializing to a 
quadratic potential V(0) = 0^). Especially the sign of this combination is important since 
a large negative value can push the bounce in the deep quantum regime even for large 
H. We have Ci — C2 = 2{AHqY > for the energy fluctuation of free solutions. If we 
use A ^ Hq which, as a consequence of the reality condition, is satisfied for free solutions 
with large Hq and assume unsqueezed states C3 ~ C4, C5 ~ Cq and put the extra condition 
C3 ~ C5, the remaining terms cancel and Gf"^ — G^ is positive. But for squeezed states 
there can be significant contributions from fluctuations at the bounce, which can even make 
G'l'^ — G^ negative there. In such a situation, it is important to know the precise state, 
i.e. all parameters ci, C2, ai, 02, and ^2, in order to determine the quantum nature of 
the bounce. However, the relevant squeezing parameters cannot be fully determined from 
using the state at only one side of the bounce [29j due to exponential suppression factors 
of some of the integration constants in fl29|) - fl32|) if one restricts to a fixed sign. Thus, 
the precise quantum nature of the bounce may always depend on what assumptions one 
makes for a state. 



5 Conclusions 

Loop quantum cosmology implements the discrete nature of quantum gravity, in contrast 
to Wheeler-DeWitt quantum cosmology. Its dynamics is determined by a difference, rather 
than differential equation, resulting in singularity-free solutions by general arguments in 
symmetric models [3, E] • But in general no intuitive geometrical picture for the avoidance 
of the classical big bang singularity is provided; typically the regime around a classical 
singularity is of a highly quantum nature. Nevertheless, such regimes may be probed by 
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effective equations, a first example being presented in [TO]. Surprisingly, a very smooth 
bounce picture resulted in the model studied there (an isotropic cosmology sourced by a 
free, massless scalar), with the whole history described well by simply replacing a connec- 
tion component c in the Hamiltonian constraint by sine. This replacement is motivated 
by the use of holonomies, rather than connection components themselves, in loop quantum 
gravity. But there are undoubtedly other corrections in a general quantum system as well 
as other types of corrections expected from a loop quantization, and it was unexpected that 
only the holonomy replacement provides effective equations in a precise manner. These 
equations were indeed shown to be effective equations in a strict sense since other quan- 
tum corrections in this specific model are simply absent. But this analysis also showed that 
the validity of effective equations obtained by the holonomy replacement is very special 
and tied to the free scalar model. 

In this paper, we have provided a detailed systematic analysis for massive or self- 
interacting scalars. The situation turns out to be remarkably different from the free scalar 
models: back-reaction of fluctuations on expectation values cannot be ignored, giving rise 
to extra quantum corrections to the classical equations. Moreover, standard as well as more 
involved choices of setting up an adiabatic approximation of fluctuations were unsuccessful 
in that the adiabaticity assumption was in conflict either with uncertainty relations or did 
not lead to any complete solution to all equations of motion. Thus, for the cases analyzed 
here there is no analog of a low-energy effective action. There may be regimes where an 
adiabaticity assumption is consistent, but this had to be rather contrived. Our conclusion 
is that in general only higher dimensional effective systems exist for quantum cosmology, 
which crucially contain independent quantum degrees of freedom. In particular, the holon- 
omy replacement of c by sine in loop quantum cosmology is not sufficient, although it still 
occurs due to the use of J = pe"^ as basic variable in the setup of effective equations. 

In general, the holonomy replacement by itself overlooks crucial quantum effects caused 
by back-reactions of a spreading and deforming state on its expectation values. (This can 
be seen as a crude way of writing effective equations by setting = 0, which obviously 
violates uncertainty relations.) This is important to realize since the replacement has 
often occurred in the recent literature, see e.g. [301 El Ell E21 ESj. It is used either in 
the form of sine in the Friedmann equation, or by using equations of motion to solve 
sine in terms of d, which gives rise to a correction term of energy density squared to 
the classical Friedmann equation [30]: starting from a Hamiltonian constraint of the form 
C = —fj.'"^ siii^ {f^c)^^/\p\ + |pp''^Pmatter (whcre n may be a ]9-dependent function) one derives 
the Hamiltonian equation of motion p = 2//~^ sin(/ie) cos(/ie)-\/[p|. From \p\ = a^, on the 
other hand, we have p = 2aa such that the constraint C = can be expressed as the 
corrected Friedmann equation 

a2 = ^-2 sin2(/ie)(l - Sin^ifxc)) = aVmatter(l - (/Ua)Vmatter) (US) 

using C = 0. This equation looks very suggestive and easily generalizable to arbitrary mat- 
ter densities, which in fact has often been done in the recent past. However, this equation 
is applicable only to an exactly isotropic model sourced by a free, massless scalar. Even 
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before embarking on a systematic analysis of effective equations as in this paper, the holon- 
omy replacement sin c could not be expressed simply as a p^-correction (or possibly higher 
power) if anisotropics or inhomogeneities are included. Then, several independent gravita- 
tional variables exist and one cannot simply solve for a single a to bring the gravitational 
part of the Friedmann equation in classical form and have only the matter contribution 
corrected. 

Moreover, our analysis has shown that the holonomy replacement is not the only quan- 
tum correction when matter interactions are allowed for. (Similar effects happen when 
anisotropics or inhomogeneities are included in addition to the complications for solving 
for a mentioned above.) Quantum fluctuation terms are not only to be included in the effec- 
tive equations but even to be kept as independent variables^ The replacement of c by sin c 
by itself does not result in proper effective equations at all; such equations are rather to be 
considered phenomenological by isolating one quantum effect but ignoring others. (A sim- 
ilar type of phenomenological equations without proper effective analysis occurred in [M] 
as the first such example in loop quantum cosmology.) Studies of such phenomenological 
equations are valuable, but to avoid misunderstanding they should be clearly designated as 
"phenomenological" rather than "effective" unless an analysis of proper quantum aspects 
has been performed. 

Compared to corrections from inverse powers of metric components, which were utilized 
in [31] and elsewhere, for higher order corrections in c as they result from the holonomy 
replacement it is more important to keep fluctuation terms. Both are related to higher 
curvature corrections (higher powers of a time derivatives of metric components on the one 
hand and higher order time derivatives on the other) as they are usually expected from 
effective actions. These corrections are better not to be separated, or else violations of 
covariance may easily occur. 

Proper effective equations for flat isotropic models sourced by massive or self-interacting 
scalars have been provided here to first order in h, provided by moments of second order. 
We have seen higher powers of c playing a role in the variable J, which also enters equations 
of motion. In addition, there are clearly quantum degrees of freedom corresponding to 
independent fluctuations in the effective equations. This is in agreement with expectations 
from higher curvature terms, although not in obvious correspondence. In our case, we 
have true independent quantum degrees of freedom even at the perturbative level, while 
higher derivatives in effective actions do not provide new degrees of freedom as independent 
solutions in a consistent perturbation scheme. 

We have analyzed the effective equations perturbatively, making sure that reality condi- 
tions are respected. The solutions correspond to moments computed from a physically nor- 
malized state. At first perturbative order, there is no strong back-reaction for a Wheeler- 
DeWitt quantization, or for most of a single collapsing or expanding branch in a loop 
quantization. But quantum back-reaction effects are noticeable for squeezed states around 

■^Even if they were included in the effective potential, matter terms themselves would be c-dependent, 
and (|113p would only be an implicit equation for d whose properties would differ from an explicit one with 
a a-independent right hand side in (|113p . 
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the bounce of loop quantum cosmology due to the influence of the growing branch on the 
collapsing one. Since squeezing parameters of a bouncing state are not determined by 
properties in one of the two large volume regime only [29], no general statement about 
the presence of a semiclassical bounce is available. There are quantum parameters whose 
influence is negligible at large volume but which will become important near a bounce and 
determine its fate. This shows the caution required for statements about the behavior 
of a bounce even in a perturbative regime around the exact bounce model. While a full 
analysis is beyond the scope of this paper, we have provided the setting based on a proper 
derivation and analysis of effective equations. 
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A Poisson brackets between fluctuations in loop quan- 
tum cosmology 

Here we provide more information on the derivation of Poisson brackets (178|) - (!86|) . All 
derivations follow the same scheme, except for ( ITHll which does not give rise to third order 
moments because [JJ\p^] = = [J^J, p^]. In this case, we have 



For the remaining brackets there are contributions by moments of third order such as 




Qp. 




where we used JpJ = \{pJ'^ + J^; 



''p). We will also use the third moment 
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obtained by a similar calculation, or just by exchanging p and J in G^'^^ . (But note that 
the final step in fllMp does differ since we now have 

pJp=^{p{pJ-hJ) + {Jp + hJ)p) = ^{p^J + Jp^-h[p,J]) = ^{p^J + Jp^-!fj). (116) 

This allows us to express the Poisson brackets between fluctuations in terms of moments 
and expectation values: 

{GPP,G'^J} = l{[p^,.p])-lj{[p^^j])-lp(ip^.p]) + {p^^j^} 

= -2i{pJ^ + Pp) + 2iJ{pJ + Jp) + Aip{J^) - AipJ^ 
= -AiG^^-^ -AiJG'P-^ 

using (11141) in the last step. Furthermore, we obtain 
{G^G^-^} = -Li^[p\pUjp])-^p{[p\J])- 

= -^'^{P^J + 2pJp + Jp^) + 2ip{pJ + Jp) - 2ipV 

2 - - - 1 

= — + pJp + Jp'^) + -ih^J + Aip{G^-^ + pJ) - 2ip^J 

3 6 

= -2iGP^^ -2iJGP^ +-ih^J 

6 

where we used (1116^ to bring the third moment in the form as it appears in (lllSp . The 
remaining brackets follow from similar calculations, although they can become more lengthy 
if all three basic operators p, J and are involved. Skipping some of the details of 
calculational steps already encountered, and occasionally using the commutation relations 
as well as J = p^ and the associated reality condition fl26|) . we further have 

{G'J^ G'J] = l-l^jp + Jt J, P]) - ^J{[JJ^ + Jt J, J]) - ^J{[J\ P]) + { JJ, J2} 

= -i{.Pp + 2 JpJ + p.P + 2h.P) + AiJ{pJ + Jp) - Aipj'^ + 2inj'^ 
= -AiG^^'^ - SiJiG^-^ + pJ) - Aip{G^'^ + JJ) + SipJ^ - 2ih{G-^-^ + J^) 
+8i J(G^-^ + pJ) - AipJ^ + 2ihJ^ = -AiGP'^-' - Ai{p + h/2)G^'' , 

{G'^G^-'} = l^{[Jp + j\j,pJ + Jp])--^p{[JJ^ + j\j,J]) 

Ain 2ifi 2in 

= -^z(]3^J + 2]3Jp + Jf) + z(p + ;i)(pj + Jp) + 2ii;ipJ + 2iJl~^f) 

+ Uj{J.P + J"^ J) - iJ{.P) - 22p(p + fl/2) J 

= -2iGPP-^ - 2z(p + hl2)G''-^ + z JG-^-^' - i JG^-^ + -ih^J , 

6 
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= 2i{J^) - AiJ{j'^) + 2iJ^ = 2iG-^-^-^ + 2iJG-^^ , 

{G'-^,G^'~} = ±{[j\pJ^ + J^p])-}^J{[J^,p])-^J{[J,pJ^ + Pp]) + {J^^ 
2tn in in 

= 2i{p^J + pjp + Jp^ + 2h{pJ + Jp)) - 2ip{pJ + Jp + hJ) - 2i J( J2) 
-iJ{J.P + J"^ J + Ap{p + h/2)) + 2i J + Aip{p + ^2) J 

= QiG^^-^ + 8i(p + n/2)GP'^ + 2iJGPP - 2iJG^^ - 2iJG^^ + 2ip^J - 2iJ^J 
+2ihpJ 

= GiGPP^ + 8i{p + h/2)GP^ - 2iJG^^ - iff J, 
{GP'^ G^'} = ^{[pJ + Jp,pJ^ + J^p]) - ^P{[lpJ^ + J^V]) - ^JilP^PJ^ + J^P]) 

"2^^^^^^^ ^^^^ " i^HiPJ + Jp,p\) + {pJ.pJ} 

= \^8p^ + P{JJ^ + 2J^ J) + {IP + 2.Pj)p + Jp.P + .Pp.! + Ahp^) 

-ip{Jp + Pj + 4]3^ + 2hp) - ^iJ{pP + Pp) - ^iJ{pJ + Jp) 

+2ipJJ + 2ip^{p + n/2) 
= AiG^PP + i{8p + 3h)GPP - 2ipG'^^ - iJG^-^' - iJG''-^ + 2ip^ - 2ipJJ 

+2ihp^ + -ih^p + \ih^ 
2 4 

= AiG'P'P'P + Qi{p + h/2)GPP - iJG'P^' - iJG^^ + ]^ih^p + hih^ 

and 

{G'^\G'-^} = Uj\{Pr])-^j{[j,{Pr])-^j{[j',p])+{.p,.p} 

in in in 

= 2i{.JpP + .Jpp + pp. J + Pp.! + h{!P + P.!)) - 4iJ{Pp + pp + hp) 

-AiJ{Jp + p! + hJ) + 8i{p + h/2) J J 
= StGPPP + 16i{p + h/2)GPP + 8i{p + h/2)G^^ - SiJG^^' 
-SiJGP-^ + Aih^p + 2ih^ . 
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